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Abstract 

An algebraic method especially suited to describe strongly anharmonic 
vibrational spectra in molecules may be an appropriate framework to study 
vibrational spectra of NaJ clusters, where nearly flat potential energy sur- 
faces and the appearance of close lying isomers have been reported. As an 
illustration we describe the model and apply it to the Be4, HjJ", Be3 and NaJ 
clusters. 



INTRODUCTION 

The study of metallic clusters has been mainly concerned with the electronic 
properties of these systems. It has been suggested, however, that in analogy to the 
way crystalline structure determines the optical response and metallic properties in 
solid-state physics, the geometric structure may play a significant role in the case 
of cluster physics jl], ||] . Small finite size systems present special problems which 
cannot directly be dealt with by means of the methods applied to the bulk. Nu- 
clear theorists, on the other hand, have much experience dealing with "the many 
body problem" and the application of these methods to clusters is already lead- 
ing to significant contributions to their understanding. These techniques, however, 
usually ignore the position of atomic nuclei H. In Jl|, ab initio calculations for 



Na+(n = 2 — 9, 11, 21) clusters at low temperature are reported, which allow the 
assignment of specific cluster geometries by comparing the theoretical results with 
experimental optical depletion spectra for these systems. While for dimers and 
trimers the rovibrational spectra lead to a precise determination of internuclear 
distances and ground-state potential energy surfaces (P.E.S.), up to now the reso- 
lution of the vibronic structure for n > 4 is not enough to carry out a theoretical 
analysis that would confirm the results of ref. [2] and give information on the P.E.S. 
for the heavier clusters wa . Even if the rovibrational data become available in the 
near future, the calculations of [2] imply that the corresponding potential surfaces 
are in general quite flat and thus one expects anharmonic effects to be significant. 
Also, in a number of cases close lying isomeric structures are found H. Since the 
usual molecular physics calculations are well suited for deep minima and not for 
"soft" molecules, it is important to consider alternative models in order to analyse 
the experimental results. In this paper we describe a new model which was designed 
to incorporate anharmonic behavior from the outset and which has been applied to 
a number of molecules, including the Naf cluster Q-Q. We illustrate the method 
with applications to some D^h molecules and clusters and to the Be4 cluster. 



THE U(2) VIBRON MODEL 

The model is based on the isomorphism of the U(2) Lie algebra and the one 
dimensional Morse oscillator, whose eigenstates can be associated with U(2) D 
SO{2) states ||. In the framework of the model the total number of bosons N 
is fixed by the potential shape and the eigenvalue m of the SO(2) generator J z , 
takes the values m = ±N/2, ±(N — 2)/2, . . .. The Morse spectrum is reproduced 
twice and consequently for these applications the m-values must be restricted to 
be positive. In terms of the U(2) algebra, the Morse Hamiltonian has the algebraic 
realization 

H = A{^-Jl) . (1) 

with eigenvalues 

E M = AN[(v + 1/2) - v 2 /N] , (2) 

where the label v = j — m denotes the number of quanta in the oscillator. The 
parameters N and A appearing in (2) are related with the usual harmonic and 
anharmonic constants uj e and x e ui e used in spectroscopy fic|| . We now consider 
the Ui(2) D SU l {2) D SO l (2) algebra, generated by the set {G t } = 
{Ni, J+.i, J-.i, Jo,i}, satisfying the commutation relations 

[Jo,i,J±,i] — ±J±,i ) [J+,i>3-,i] — 2 Jo,i , [Ni,J^i] — 0, (3) 

with [a = ±,0. For the symmetric irreducible representation [iVj,0] of Ui(2) one 
can show that the Casimir operator is given by || 3 2 — Ni(Ni + 2)/4, from which 
follows the identification j$ — Ni/2. 



In the algebraic approach each relevant interatomic interaction is associated 
with a Ui(2) algebra pi ||. As a first example, we consider the Be4 cluster, which 
has a tetrahedral shape. T> 3 h molecules can be similarly treated. In the Be4 case 
there are six Ui(2) algebras involved (i = 1, ... ,6). The operators in the model 
are expressed in terms of the generators of these algebras, and the symmetry re- 
quirements of the tetrahedral group Ti can be readily imposed || O]. The local 
operators {G,} acting on bond i can be projected to any of the fundamental irreps 
r = A\, E and F 2 . Using the generators (3) we obtain the Ti tensors 

6 
z=l 

where [i = ±, and 7 denotes the component of T. The expansion coefficients are 
the same as those given in the one phonon wave functions 0. The Hamiltonian 
operator can be constructed by repeated couplings of these tensors to a total sym- 
metry Ax, since it must commute with all operations in Td. This is accomplished 
by means of the 7rf-Clebsch-Gordan coefficients || O, |TJ| . 

All calculations can be carried out in a symmetry-adapted basis, which is pro- 
jected from the local basis 

Ux(2) <g> ••• ® U 6 (2) D SOx{2) ® ••• ® S0 6 (2) D 50(2) 

I II I I (5) 

I [Nx] , ... , [JV 6 ] ; vx , ... , v 6 ; V ) 

in which each anharmonic oscillator is well defined. By symmetry considerations, 
Ni — N for the six oscillators, «j, denotes the number of quanta in bond i and 



V = J2i v i i s the total number of quanta |13, g]. The local basis states for each 
oscillator are usually written as \Ni, Vi), where v, L = (Ni — 2mj)/2 = 0, 1, . . . [iVj/2] 
denotes the number of oscillator quanta in the i-th oscillator. The phonon states 
\ v <p^ > can be constructed using the Clebsch-Gordan coefficients of Td J|, Q. Since 
all operators are expressed in terms of powers of the E/»(2) generators, their matrix 
elements can be easily evaluated in closed form. The symmetry-adapted operators 
(11) and states are the building blocks of the model. 

We now proceed to explicitly construct the Be4 Hamiltonian. For interactions 
that are at most quadratic in the generators the procedure yields 

H = wi Ha x + u 2 He + ^3 Hf 2 + b 2 Ve + b 3 Vf 2 , (6) 

with 

7 

Vr - ^E^ 7 ^ 7I (7) 



N ■ 
7 



The five interaction terms in Eq. (10) correspond to linear combinations of the ones 
obtained in lowest order in |6[ However, it is necessary to include interactions 
which are related to the vibrational angular momenta associated with the degen- 
erate modes E and F2 ■ These kind of terms is absent in the former versions of the 
model ||, We now proceed to show how they can be obtained in the present 
model. In configuration space the vibrational angular momentum operator for the 
E mode is given by Q 

where qf and qf are the normal coordinates associated to the E mode. This rela- 
tion can be transformed to the algebraic space by means of the harmonic oscillator 
operators 



to obtain 



b»-b^b?) . (10) 



Here = i bi, with a similar form for b^ while the 

(4). In order to find the algebraic expression for l A2 we first introduce a scale 
transformation in (11) 

b\ = J-.i/v 7 ^, bi = J+,i/VNi. (11) 

The relevant commutator can be expressed as 

Ml] = ±-[J+,i,J-d = ^2J , 4 = 1-^ , (12) 

where 

- ^ f fi»\ 

Vi = — - Jo,i ■ (13) 

The other two commutators in (11) are not modified by (11). In the harmonic limit, 
which is defined by Ni — ► 00, Eq. (12) reduces to the standard boson commutator 
\bi,bl) = 1. This limit corresponds to a contraction of SU(2) to the Weyl algebra 
and can be used to obtain a geometric interpretation of the algebraic operators in 
terms of those in configuration space. In the opposite sense, Eq. (5) provides a pro- 
cedure to construct the anharmonic representation of harmonic operators through 
the correspondence b\ — > b\ — J-^/y/Nl and bi — > bi = J+^/\fNl. Applying this 
method to the vibrational angular momentum (4) we find 

i M - -(i-Vf.-lVfJ • (14) 



Table 1: Vibrational excitations of Be4 using the algebraic Hamiltonian with 
parameters given in the text. The ab initio (N — ► oo) spectrum is generated with 
the parameters from [ jT5[ . The energies are given in cm -1 . 
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For the vibrational angular momentum l^ 1 associated with the F 2 mode we find a 
similar expression. 

We can now use our model to fit the spectroscopic data of several polyatomic 
molecules. In the case of Be4 the energy spectrum was analyzed by ab initio meth- 
ods in Jl5| , where force- field constants corresponding to an expansion of the po- 
tential up to fourth order in the normal coordinates and momenta were evaluated. 
We have generated the ab initio spectrum up to three phonons using the analysis 
in juj . For the algebraic Hamiltonian we take Q 

H = lu\ Hai + ^2 'He + W3 Hf 2 + X 33 (Hf 2 ^J + -X12 T^Al J 

+X 13 (h Ai Hf 2 ) + 333 Z T + *33 033 + t 23 023 ■ (15) 

7 

The terms O33 and O23 represent the algebraic form of the corresponding interac- 
tions in juj which are responsible for the splitting of the vibrational levels in the 
(1/1,^,14) = (0,0°, 2 2 ) and the (0,! 1 ,! 1 ) overtones. 



Table 2: Least-square energy fit for the vibrational excitations of H3 , Be3 and Najj 
The energy differences AE = E t h — E exp are given in cm 
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Note that the Be4 Hamiltonian (15) preserves the total number of quanta V. 
This is a good approximation for this case according to the analysis of |l4|, 
but it is known that Fermi resonances can occur for certain molecules when the 
fundamental mode frequencies are such that (V, V) states with V 7^ V are close 
in energy. These interactions can be introduced in the Hamiltonian by means of 
a polyad analysis fl6|| . For T> 3h molecules we can follow an analogous procedure, 
namely, we can construct the T)^ symmetry-adapted operators and states corre- 
sponding to (4) and and carry out the building up procedure to construct the 
Hamiltonian and higher phonon states, using in this case the appropriate projection 
operators and Clebsch-Gordan coefficients M, Ww. 



EXAMPLES 



We now present the results of our least-square fits to the energy spectra of Be4, 
Be 3 , NaJ and H^". In Table 1 we show the fit to Be4 using the Hamiltonian (15). 
The fit includes all levels up to V = 4 quanta and gives a r.m.s. deviation of 2.6 



cm -1 , which can be considered of spectroscopic quality. In Table 1 we only show 
the results for the levels with V < 3. We point out that in |l4|, [l5| several higher 
order interactions are present which we have neglected. Since our model can be 
put into a one to one correspondence with the configuration space calculations, it 
is in fact possible to improve the accuracy of the fit considerably, but we have used 
a simpler Hamiltonian than the one of |l4|, [l5|]. When no ab initio calculations 
are available (or feasible) the present approach can be used empirically, achieving 
increasingly good fits by the inclusion of higher order interactions . In Table 2 
we present fits to the spectra of Be3, Najj" and H^" up to three phonons. While 
remarkably accurate descriptions of the first two molecules can be achieved using 
a four-parameter Hamiltonian, we were forced to include four additional higher 
order terms in the H^" Hamiltonian in order to properly describe this molecule. 
This is in accordance with the work of Carter and Meyer Jl7| , who were forced to 
include twice as many terms in the potential energy surface for H^" than for the 
Na^ molecule. The H^" ion is a very "soft" molecule which, due to the light mass of 
its atomic constituents carries out large amplitude oscillations from its equilibrium 
position s |l7[ . This may also be the case for the potentials associated to metallic 
clusters fell. 



SUMMARY 



In this paper we have studied the vibrational excitations of several molecules 
and clusters in a symmetry-adapted algebraic model. These studies suggest that the 
symmetry-adapted algebraic model provides a numerically efficient tool to study 
molecular vibrations with high precision. The main difference with other methods is 
the use of symmetry-adapted tensors in the construction of the Hamiltonian. In this 
approach, the anharmonicity can be introduced from the outset, interactions can be 
constructed in a systematic way, each term has a direct physical interpretation, and 
spurious modes can be eliminated exactly [ fl6[ . It is important to further explore 
the scope and applicability of the present approach, and in particular study the 
vibrational spectra in metallic clusters, including the appearance of isomeric shapes 
with different point symmetries. 
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